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ABSTRACT 


^A  general  solution  of  the  heat  diffusion  equation  is  presented  for  the 
case  of  a  semi-infinite  two-^l^ye^  system  that  is  heated  with  a  localized  cw 
laser  beam  of  Gaussian  or  uniform  circular  intensity  profile.  As  an  example, 
this  theory  is  applied  to  thin  layers  of  a  thermoluminescent  material  on  a 
substrate.  The  transient  thermoluminescence  emission  response  is  calculated 
and  compared  with  experimental  results.  Illustrating  the  validity  of  the 
solutions.  Applications  to  transient  laser  heating  of  oxide/semiconductor 
sandwiches  are  discussed. 

The  theory  of  heating  a  semi-infinite  two-layer  system  with  a  localized 

P. 

CO^  laser  beam  of  uniform  circular  intensity  profile  is  also  applied  to 
configurations  consisting  of  thin  thermoluminescent  LlF;Ti,Mg  phosphor  layers 
on  borosilicate  glass  substrates. ■  We  study  thermoluminescence  because  it 


serves  as  a  convenient  monitor  the  temperature  distribution  and  because 
of  its  Importance  in  developing  the  laser  heating  technique  for  solid-state 
dosimetry.  Experimental' and  calculated  results  are  compared  in  an  attempt  to 
completely  characterize  the  thermoluminescence  response  curves.  The  two-laye 
system  is  heated  in  two  different  modes:  the  laser  beam  impinges  onto  (a) 
the  phosphor  layer,  and  (b)  the  glass  substrate.  We  investigate  in  detail 
changes  in  the  thermo luminescence  response  due  to  different  laser  powers  and 
variations  in  the  optical  and  thermal  properties  of  the  samples. 

J) 

Two  optically  stlmulable  phosphors  (MgS:Ce,Sm  and  CaS:Ce,Sm)  have  been 
developed  with  the  goal  to  eventually  apply  them  to  imaging  of  neutron  dose 
distributions.  These  phosphors  exhibit  luminescence  emission  spectra  that 
are  different  for  UV  excitation  and  for  IR  stimulation,  a  property  necessary 
for  nearly  real-time  radiography. 
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RESEARCH  GOALS 


The  purpose  of  chis  research  effort  is  the  investigation  of  the  basic 
characteristics  of  two  differnt  types  of  laser-stimulated  luminescence 
phenomena  in  view  of  potential  applications  in  the  dosimetry  of  ionizing 
radiation  and  of  radiation  imaging. 

The  first  of  these,  laser-heated  thermoluminescence,  was  the  main  subject  of 
research  in  the  contract  period  from  7-15-83  through  3-15-86.  A  study  of  the 
potential  use  of  optically  stimulated  luminescence  for  real-time  imaging  of 
neutron  dose  distributions  has  been  phased  in  since  5-25-86. 

The  heat  transfer  mechnanism  from  a  CO  laser  beam,  having  a  circular  uniform 
intensity  profile,  to  a  semi-infinite  thin  thermoluminescence  phosphor  layer 
and  the  resulting  luminescence  response  has  been  studied  in  detail. 
Experimental  and  computational  methods  were  developed  for  the  chacterization 
of  a  given  thermoluminescence  dosimeter  configuration.  All  its  relevant 
optical  and  thermal  properties  can  be  determined  by  these  methods.  Based  on 
this  knowledge,  the  expected  thermoluminescence  response  can  be 
computer-simulated  as  a  function  of  a  number  of  experimental  parameters.  This 
in  turn  provides  a  tool  to  optimize  the  dosimeter  design  for  applications  in 
dosimetry. 

Initial  experiments  are  discussed  which  we  performed  with  a  group  of  Ce/Sm 
co-doped  magnesium  and  calcium  sulphides  that  can  be  stimulated  with  1.064  urn 
photons  from  a  Nd;YAG  laser.  The  goal  is  to  develop  phosphors  that  exhibit 
distinctly  differnt  luminescence  emission  spectra  during  excitation  with 
ionizing  radiation  as  compared  to  those  emitted  during  infra-red  stimulation. 
This  characteristic  property  is  essential  for  real-time  or  nearly  real-time 
radiographic  imaging. 
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ABSTRACT 

A  general  solution  of  the  heat  diffusion  equation  is  presented  fnr  the  case  of  a 
senii-inilnite  two-layer  system  that  is  heated  with  a  localized  cw  laser  beam  of  Gaussian  or 
uniform  circular  intensity  profile.  As  an  exampte,  this  theory  is  applied  to  thin  layers  of  a 
thermoluminescent  material  on  a  substrate.  The  transient  thermoluminescence  emission 
response  is  calculated  and  con:q)ared  with  experimental  results,  illustrating  the  validity  of 
the  solutions.  Applications  to  transient  laser  heating  of  oxide-semiconductor  sandwiches 
are  discussed. 


INTRODUCTION 


The  interest  in  heating  senii'infinite  slabs  with  localized  laser  beams  stems  from  the 
j^lication  of  laser  heating  in  tbeimoluminescence  dosimesy ,  which  recently  has  received 
considerable  attention  [1-5].  The  eariy  week  in  tiiis  field  was  restricted  to  the  use  of  cw 
CO2 beams  having  Gaussian  intensi^  profiles.  The  tiieoretical  approach  talcen  for 

the  calculations  of  the  thermoluminescence  (TL)  response  of  a  thin  TL  layer  oh  a  glass 
substrate  assumed  that  the  TL  layer  thickness  is  negligibly  small  compared  to  that  of  the 
substrate.  All  thermal  properties  of  the  sample  were  assumed  to  be  due  to  the  substrate 
alone.  Only  the  reflectivity  and  absoptance,  due  to  modification  of  the  glass  substrate  by 
the  thin  TL  layer,  were  taken  into  account  The  agreement  between  experimental  TL 
responses  and  calculated  ones  was  satisfactory  for  a  number  of  layer  configurations 
consisting  of  LiF:Ti,Mg  powder  (Harshaw/Filtrol  TLD-100)  on  borosilicate  glass  [5]. 

For  practical  applications  of  laser  stimulated  thermoluminescence  in  dosimetry  it 
became  apparent  however,  that  Gaussian  beam  profiles  impose  unnecessary  restrictions 
[6].  We  have,  therefore,  developed  the  optical  technology  to  produce  spatially  uniform 
beam  intensity  profiles  [7]  which  was  modified  for  uniform  circular  beams,  and  we 
subsequently  attempted  to  calculate  the  thermoluminescence  response  functions  obtained 
with  such  a  beam  on  the  basis  of  the  simple  one-layer  tbemy  {xesented  in  Ref.  5.  It  turned 
out  that  this  theory  did  not  yield  the  characteristic  TL  glow  curves  measured  under  these 
conditions.  This  failure  prompted  us  to  work  out  the  theory  of  localized  laser  beam  heating 
of  a  two-layer  system.  It  is  the  purpose  of  the  present  paper  to  present  the  general 
solutions  cf  the  heat  diffusion  equation  for  this  case  and  apply  the  theory  to  the  same  Qpe 
of  thin-layer  thermolummescence  dosimetry  configuration  (LiF;Ti>Ig  TL  powder  on  a 
borosilicate  glass  substrate)  used  for  the  work  described  in  Ref.  5.  For  completeness  we 
have  also  recalculated  the  transient  heating  and  the  resulting  TL  response  of  this  two-Iaya 
configuration  for  Gaussian  heating  beams.  It  is  perhaps  worthwhile  to  point  out  that  the 


7 

solutions  described  in  this  paper  are  of  general  validity  and  may  well  be  useful  in  other 
fields  of  endeavor;  Le.,  laser  annealing  of  semiconductors. 

n.  SOLUTIONS  TO  THE  HEAT  EQUATIONS 

The  system  we  consider  consists  of  two  layers  of  material  in  perfect  thermal 
contact  As  an  example  picture  again  a  thin  layer  of  a  thermoluminescent  phosphor  on  a 
glass  substrate  or  an  oxide  layer  on  a  semiconductor.  Each  layer  is  characterized  by  its 
abstxption  coefficient  for  the  laser  photons,  thermal  conductivity  Iq,  thermal 
diffusivity,  ccj,  and  thickness,  Ij.  The  index  i «  1,2  denotes  the  respective  layer  and  the 
total  thickness  is  L  ^  ij  + 12.  We  also  assume  that  the  laser  beam  impinges  along  the  z- 
axis  perpendicularly  to  the  x,y  plane  of  the  two-layer  configuration.  The  laser  beam  profile 
is  either  circular  and  uniform  (Fig.  1)  or  Gaussian. 

The  thermal  diffusion  equation  in  cylindrical  coordinates  (ignoring  the  effect  of  the 
finite  speed  of  heat  propagation  [8]}  is  given  by 

rl3(rki  dTi(r,z,t)/9r]/8r  +  32T(r,z,t)/3z2  +  gj(r,2,t)  «  (ki/ai)dTi(r,2,t)/9t  (1) 

We  solve  it  for  the  following  boundary  conditions: 

1 .  the  two-layer  slab  is  infinite  in  radial  direction  and  T(r,z,t)  -^0  as  r 

2.  no  heat  loss  occurs  at  the  boundaries;  i.e.,  dTj/dz  »  0  at  z  *  0  (i » 1)  and 
z-L(i-2); 

3.  at  the  interface  (z  » Ij),  Tj  «  T2  and  kj  ST^/dz*  k2 

The  initial  condition  is  Ti(r,z,t) «  Tg;  Le.,  room  temperature,  when  the  laser  is 
turned  on  at  t «  0.  The  temperature  AT|(r,z,t) »  T|(r,z,t)  -  Tg,  obtained  from  Eq.  1,  is  the 
actual  temperature  increase  above  Tg. 

Correcting  for  reflection  loss  at  the  front  surface  at  z  >■  0,  the  source  functions 
become  for  0  t  ^  tp  (q>  is  the  laser  pulse  length): 


gi(rAt)  »  Pi(l  -  Ri)(P*/7to)2)exp(-piz)f(r,o)) 


(2) 
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and 

g2(rAt)  -  Mjd  -  Ri)(l  -  R2)  (P*/JC<o2)exp[-Ui/i  -  (3) 

where  f(r,a)) «  H(a>>r)  for  the  unifonn  beam  and  f(r,(ii>)  *  exp(-i2/cii>^  for  the  Gaussian 

beam.  P*  is  the  total  laser  power  (not  the  power  density)  incident  on  the  front  surface,  H 

o 

the  Heaviside  step  function,  R>  the  reflectivity  layer  i,  and^the  radius  of  the  uniform 
circular  beam  or  the  1/e  half-width  for  the  Gaussian  beam. 

We  solve  Eq.  (1)  with  the  initial  and  boundary  conditions  given  above  by  the 
Green's  function  technique  [9].  The  general  solution,  written  compactly,  is  then: 

ATi(r,z,t)  «  !q  dtj^  r'dr'[(ai/ki)  jJ*  d2'Gii(r,2,t/r’,z’,t)g,(r’,2',x)  + 


(Oj/kj)  dz’Gi2(r,2,t/r’,2’,T)g2(r’,2’,T)l  (4) 

Gjj  »  (kj/Oj)  Q  (Y,<lYi/otj)«xp(-(Yi2  +  ni^(t*T)]Jo(Yir/^ai)  JoCY/Z^Oj) 

[Xi(2)Xj(2')]/N  (5) 

2 

N  «  ^!^ki/i/2cxi(l  +  sin29i/29i]/cos2%)  where  9i»§i/i/V(Xj  (6) 

/^cos(5i2/Vai)/cos(^i/i/Vai)  for  i^l 

Xi(2)  -  (7) 

(cos(42(L-z)/Va2)/cos(§2V“^ot2)  for  i»2 

From  the  bourxlaiy  conditions  at  the  interface  we  have 

Yi/Voi  »  s  p. ,  (8) 

i  kj^i/Votj  tanC^i/j/Voi)  -  0  (9) 

1-1 


and 


W  *  §1^  +  p^(ai  -  “2)  • 


(10) 
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The  rise  in  temperature  dien  becomes 

ATi(rAt)  -  [(l-Ri)P*/)ca)2N]  dt  r’dr’  pdp  exp[aij32(t.x)]Jo(pr)  x 
Jo(pr')f(r’.a))  Xj(2)  exp(-^2(t-t)]{p.j2  exp(-|Xi/i) 

[exp(4i/i)  +  (§^/^tl^/al)sin(§^/lWal)  - 

COS(§i/i/Vai)]/[(^Li2+4i2/otj) 

cos(§i/i/Vaj)]  +  (l-R2)/p.2^ 

CXp(-4j/,)[C0S(§2V^Ct2)  +  (42W«)S»n(§2V^“2)  ' 

CXTp(r\S^lT)V[{\l-^  +  ^22/a2)COS(^2V^«2)l}- 

This  general  solution  can  be  simplified  for  all  those  two-layer  structures  for  which 
one  of  the  layers  is  much  thinner  than  the  other  one.  As  we  will  show  below  in  several 
examples,  the  spatio-temperal  temperature  distribution  of  such  a  configuration  can  be 
calculated  by  assuming  =  02  =  (X.  This  assumption  entails  that  the  general  form  of  the 

solution  is  maintained  in  regards  to  the  deposition  of  the  laser  photons  according  to  the 
absorption  coefficients  and  II2,  however,  that  the  thermal  diffusion  is  governed  mainly 

by  the  thicker  one  of  the  two  layers.  Lateral  heat  diffusion  in  the  thinner  of  the  two  layers 
is  negligible.  Instead,  lateral  heating  in  this  layer  occurs  via  diffusion  in  the  thicker  layer 
with  heat  transport  across  the  interface  into  the  thinner  layer.  For  this  case  the  integration 
over  P  and  r'  in  Eq.  (1 1)  is  readily  carried  out,  yielding 
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,t 

ATi(r.z,t)  =  (BP-/Q)2)  Jo  dt  Q(t.t){l  +  [(1  -  Ri)/jcB] 

go 

I  cxp(-4i/i)[exp(-jti/i)cos(Tini/i)  +  (Tln|/4i)sin(Tlm/i)]/ 

n-l 

[Nn(l  Tla2/4i2)cos(Tln^^l/i)]  +  (1  -  R2) 
cxp(-Ui/i)[cos(lln,/2)  +  (Tlni/JA2)sin(Tln,/2)  •  CxpC-p^/J/ 

INn,(l  +nm^/^i2^)cos(Ttn,/2)]}Xi(z)exp[-anro^(t<)]  <12) 

with  Q(t) »  P((ii/V2at,  rHlat)  for  uniform  circular  beam  or  exp(r^/(o)2  +  4at))/ 

at 

[1  +  (4at/Q)2)]  for  Gaussian  beam  proHle,  where  P(<a^) »  Jq  o>'do>'  Io(toV)  x 
exp(-(ca'2+r2)/2]  is  called  the  circular  coverage  function  [10]  or  the  "p"  function  [1 1], 
§i/Vai  =  =  Ti„,  B  «  (1  -  Ri)[l  -  (1  -  R2)exp(-pi/i  - 112^2)-  ^2  exp(-iii/i)]/5cNo 

and,  from  the  boundary  conditions,  the  transcendental  equation 

k,  tan  tiniA  + 1^2  tan  ti„/2  »  0-  (13) 

where  the  roots  m  »  1,2, ...  are  derived  in  Appendix  A. 

We  illustrate  the  validity  of  these  results  obtained  for  both  the  uniform  circular  and 
the  Gaussian  beam  by  monitoring  the  thermoluminescence  response  of  a  thin  TL  layer  on  a 
glass  substrate.  We  choose  as  an  example  a  35-40  pm  thick  LiF:Ti>Ig  layer  on  a  150  |im 
thick  borosilicate  microscope  cover  slide.  The  transient  temperature  distribution  is 
indirectly  monitored  via  the  transient  thermoluminescence  response  of  the  thin  active  layer 
to  localized  laser  heating,  measured  with  a  photomultiplier.  We  point  out,  however,  that  in 
this  way  the  spatially  integrated  TL  signal  is  measured  without  spatial  resolution.  Only  the 
temporal  evolution  of  the  thermoluminescence  signal  is  monitored.  The  procedure  for 
calculating  the  thermoluminescence  transient  is  based  on  the  well-known  electron  kinetics 
for  this  thermoluminescent  phosphor  (see  below). 
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The  absorption  coefficient  of  LiF  for  10.6  |im  CO2  laser  photons  is  40  cm*^  and 
that  of  glass  is  very  high.  Consequently,  we  have  «  1  for  the  LiF  layer  and 

»  1  for  the  glass  substrate.  We  have  locally  heated  this  structure  with  laser  beams  of 
Gaussian  and  uniform  circular  profiles,  constant  power  P*,  and  various  beam  sizes  from 
both  sides,  that  is  the  active  TL  layer  can  either  beontopat  z»0  oratthe  bottom  at  z  » 

L,  and  we  have  monitored  the  transient  theimoluminescence  response  in  both  cases.  The 
transient  temperature  distiibudons  are  readily  obtained  from  Eq.  (12)  as  follows: 

a.  Case  1 :  thermoluminescent  layer  on  top  of  the  glass  substrate 

(y-i/i «  1) 

OD 

ATi(rAt)  -  (BP*/(D2)i5  dxQ(r,t-T){l+No  I,  cosCti^z) 

«p[-anm^  (t-t)]/[N„(l  +  r\^i\i^)  cos(Tlm/i)]}  (14) 

b.  Case  2:  thermoluminescent  layer  at  the  bottom  of  the  glass  substrate 

1) 

09 

AT2(rAt)  =  (BP’*/(o2)J5  dtQ(r,t-x){l+No  I  cos(7l„(L-z)) 

ni“l 

exp(-<xn„2  (t-x)]/[N„(l  +  Tlm2/4j2)  cos(ilm/j)  cos(ti„/2)]}  (15) 

a*  a« 

In  Eqs.  (14)  and  (15)  we  have  written  Nq  « (1/a)  Z  Iq/i  and  N,q  »  (l/2a)  Z 

C»l 

kj/jEl  +  (sin26i)/(26i)]/cos^6i,  where  9{ «  The  calculations  were  performed  with  the 

following  values  for  the  thermal  conductivities  of  the  two  layers: 

1.  TL  layer  kj  «  0.0275  W/Kcm 

2.  Glass  substrate:  k2  ■  O.CX)942  W/Kcm. 

Both  of  these  values  are  temperature  dependent  The  temperature  dependence  of  k2  for 
glass  is  given  by  k2  »  0.00942  (1  +  oATj)  with  o  «  0.00222  K*^  [12].  This  can  be 
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accounted  far  in  our  calculations  by  the  method  of  the  Kiichhoff  transformation  [5].  A 
"corrected  temperature"  TNi»To-H  <T  “‘[(1  +  2aATi)i'2  - 1]  is  introduced  with  the  same 
o.  The  temperature  dependence  of  kj  for  LiF  does  not  have  to  be  accounted  for.  This  is 
so  because  of  the  fact  that  it  is  only  a  very  thin  layer  and,  as  we  explained  above,  lateral 
heat  diffusion  is  dominated  by  the  thicker  glass  layer. 

The  thermal  diffusiviQr,  a,  assumed  to  be  identical  for  both  layers  (see  above) 
requires  special  attention.  According  m  the  glass  manufacturer  [12]  it  is  slightly 
temperature  dependent  as  well.  It  varies  from  0.0052  to  0.0057  cm^/s  between  O'C  and 
100*C  respectively,  and  it  is  expected  to  be  somewhat  higher  yet  at  the  highest  temperature 
of  4(X)*C  reached  in  our  measurements.  We  choose  not  to  correct  for  this  small  variation  in 
a  with  T.  Instead,  we  determine  an  average  value  for  the  thermal  diffusivity  by  the  best 
fit  of  the  calculated  and  the  measured  thermoluminescehce  response  curves  obtained  with 
the  unifrmn  circular  beam.  It  turns  out  that  the  high  temperature  decay  of  this  curve  is  very 
sensitive  to  slight  variations  in  a.  The  best  fit  in  all  cases  reported  on  below  (including  the 
ones  where  a  Gaussian  beam  was  employed)  was  obtained  with  a  »  0.0064  cm^/s.  This 
value  is  consistent  with  expectation  as  the  peak  of  the  TL  response  curve  occurs  at  about 
250*C,  which  appears  to  be  a  good  "average"  temperature  for  which  to  select  a  value  for  the 
thermal  diffusivity. 

For  computational  efficiency,  the  numerical  calculations  of  Eqs.  (14)  and  (15)  start 
at  the  lower  time  limit  of  t »  10^  s  rather  than  zero.  This  can  be  done  since  for  laser 
powers  that  are  believed  to  be  practical  in  thermoluminescence  dosimetry  applications  the 
temperature  rise  anywhere  in  the  sample  during  the  initial  microsecond  is  always  much 
smaller  than  the  temperature  at  which  any  thermoluminescence  signal  can  be  detected. 

The  prefactor  B  « (1  -  Ri)(l  -  R2)/JcNo  in  Eqs.  (14)  and  (15)  is  measured  by  the 

procedure  described  in  detail  in  Ref.  5.  This  procedure  determines  the  absorptance  and  B 
is,  except  for  a  constant  factor,  the  fraction  of  the  laser  power  absorbed  by  the  two-layer 
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system.  Its  measurement  automatically  includes  the  reflection  on  the  front  surface  and  the 
interface.  The  roots  obtained  from  the  boundary  condition  are  derived  in  Appendix  A. 

The  tbermoluminescence  glow  curves  we  consider  are  conunonly  called  peaks  3, 

4,  and  5.  They  are  sufficiently  well  described  by  simple  first-order  kinetics  with  attempt- 
to-escape  frequencies  of  4  x  10^^  7.3  x  10^^  s‘^,  and  4  x  lO^^  s*^  and  trap  depth  of 

1.23  eV,  1.54  eV,  and  2.17  eV,  respectively  [13].  Their  relative  intensities  are  adjusted  to 
yield  the  actual  glow  curve  measured  with  conventional  unifrnm  heating.  Peaks  1  and  2  are 
removed  by  a  10  minute  pre-anneal  treatment  at  100*C.  The  two-layer  configuration  is 
irradiated  with  ^o  gamma  rays  to  fill  the  electron  traps  and  exposed  to  a  uniform  circular 
CO2  laser  beam  of  fixed  6  J  watts  power  both  from  the  front  (active  TL  layer  exposed 

directly)  and  from  the  back  (active  layer  heamd  via  diffusion  through  the  glass  substrate, 
which  is  exposed  to  the  laser  beam).  Laser  heating  generates  a  time-dependent  transient 
spatial  temperature  distribution  which  results  in  the  TL  emission  response  shown  in  Figs. 
2(a)  and  2(b),  measured  without  spatial  resolutitm  by  a  photomultiplier  detector.  As  stated 
above,  the  one-layer  theory  does  not  yield  satisfactory  agreement  with  measured 
thermoluminescence  response  curves  when  the  heating  beam  has  a  uniform  circular 
intensity  profile.  These  curves  have  a  pronounced  peak.  However,  due  to  lateral 
diffusion,  the  thermoluminescence  intensity  does  not  vanish  completely,  when  heating  is 
continued  indefinitely  after  peak  occuxreiKe.  We  have  compared  the  simple  two-layer 
theory  and  the  one  in  which  the  temperature  dependence  of  the  thermal  conductiviQr  is  taken 
into  account  (via  the  Kirchhoff  transformation  techiuque)  with  experimental  results 
obtained  for  a  fixed  laser  power  of  watts  and  a  beam  radius  of  m  a  0.1 5  cm 
(Figs.  2(a)  and  2(b)).  The  thermal  diffusivity  chosen  was  again  a  >  0.0064  cm^/s  (an 
"average  temperature  value").  The  prefacttx:  B  was  selected  to  yield  the  best  fit  It  is  to  be 
compared  with  experimental  values,  determined  with  the  method  of  Ref.  5,  for  both  front 
and  back  heating  of  B  «  10.7  ±  25%  K-cm^/J  and  B  *  6.3  ±  20%  K-cm^/J,  respectively. 
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The  agreement  between  the  measured  and  the  calculated  thermoluminescence  response 
curves  is  very  good,  confirming  the  validity  of  the  theoretical  approach  presented  above. 

We  have  carried  out  a  detailed  program  exploring  the  thermoluminescence 
response  ctirves  obtained  under  various  experimental  conditions  for  die  purpose  of 
characterizing  this  type  of  laser  heating  in  dosimetry  applications.  The  results  will  be 
presented  in  a  forthcoming  publication  [14]. 

Figures  3(a)  and  3(b)  are  families  of  curves  measured  for  different  beam  diameters 
(1/e-widths  (O)  of  the  Gaussian  laser  beam.  The  previous  one-layer  approximation  [5] 
yielded  a  satisfactory  fit  only  by  adjusting  the  prefactor  B.  With  the  two-layer  heat 
equation  we  obtain  the  curves  drawn  through  tte  experimental  points  in  Rgs.  3(a)  and  3(b) 
in  one  case  without  correcting  for  the  temperature  dependence  of  the  thermal  conductivity 
and,  for  comparison,  when  the  temperature  dependence  is  taken  into  account  by  the  method 
of  the  Kirchhoff  transformation  as  described  above.  Again,  the  prefactor  B  had  to  be 
adjusted  slightly  for  the  best  fit  We  note  that  the  two-layer  theory  represents  a  rather 
marked  improvement  of  the  previous  <xie-layer  approximation. 

hi  conclusion,  we  have  presented  a  general  theory  of  heating  two-layer  systems 
with  localized  laser  beams  of  Gaussian  and  uniform  circular  intensity  profiles  and  checked 
its  validity  by  comparing  calculated  and  measured  thermoluminescence  response  curves 
obtained  with  thin-layer  TL  phosphor  structures  on  glass  in  a  simpler  substrate.  The 
theory  is  well  suited  for  all  systems  that  consist  of  a  thin  layer  of  a  given  material  in  perfect 
contact  with  a  thicker  substrate  of  different  optical  and  thermal  properties.  In  particular,  we 
expect  the  theory  to  be  ^licable  to  local  heating  (laser  annealing)  of  thin  oxide  or  epitaxial 
layers  on  semiconductor  or  insulator  substrates. 
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APPENDIX  A 

In  this  woric  we  have  a  pooriy  abstxhing  phosphra*  coating  »  0.2L  on  a  highly 
absorbing  glass  I2  •  0.8L  where  the  total  thickness  L  »  0.01875  cm.  In  this  case  m 
distinct  roots  'Hqi  can  be  obtained  from  the  transcendental  equation  k^  tan  6  + 
kj  tan  40  -  0  with  0  -  0.2Tin4-.  Widi  F  -  5ji/L,  h  -  3+2/k,  k  «  ki/kj,  D+  «  {h  +  [h2  - 
(1  +  4/k)]^>i^  and  D* «  {h  -  [h  +  (1  4/k)l^}*'2,  then  the  m  roots  greater  than 
tIq^O  are  given  by 

•  [Jf*  tan'‘(D*)  +  (m  -  1)1  •  F  for  m  ■  1,  6,  11,... 

[JT*  tan'^fD-^)  +  (m  - 1)1  •  F  for  m  «  2,  7,  12,... 

[1  -  JT*  tan*‘(D*)  +  (m  - 1)1  •  F  fw  m  «  3,  8,  13,... 

(1  -  JT*  tan*>(D+)  +  (m  - 1)1  •  F  for  m  »  4,  9,  14,... 

m-  F  for  m  »  5,  10,  15,... 

The  root  tIq  «  0  yields  the  first  term,  unity,  in  the  brackets  of  Eqs.  (14)  and  (15). 
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FIGURE  CAPTIONS 

Fig.  1.  Schematic  representation  of  a  unifonn  circular  beam  profile  impinging  on  the 
surface  at  z  »  0  of  a  two  layer  system  of  total  thickness  L. 

Fig.  2.  (a)TLresponsecurvemeasuredwidia6J  WCO2  laser  beam  (oa- 0.1 5  cm 

for  the  unifonn  circular  beam  profile)  by  heating  the  phosphor  coating  on  top  of  a 
glass  substrate.  The  dashed  curve  is  obtained  with  B  >  7.3  K-cm^/J,  a  »  .0064 
cm^/s,  using  Eq.  (14)  of  the  text  The  full  curve  is  the  fit  obtained  with  B  3  9. 1 
K-cm^/J,  the  same  a,  and  the  Kirchhoff  transformation. 

(b)  TL  response  when  the  laser  beam  impinges  onto  the  side  opposite  to  the 
phosphor  coating.  The  dashed  curve  is  the  fit  obtained  with  B  «  5.8  K-cm^/J  at 
W  laser  power,  a  *  .(X)64  cm^/s,  using  Eq.  (15)  of  the  text  The  solid  line  is 
the  fit  obtained  with  B  »  9.7  K-cm^/J,  the  same  oe,  and  the  Kirchhoff 
transformation. 

Fig.  3.  (a)  TL  response  curves  from  Fig.  12  Ref.  5  recalculated  using  a  4.45  W  CC)2 

laser  beam  with  a  Gaussian  intensity  ;xoEle  by  heating  the  phosphor-coated  side 
on  a  glass  substrate.  The  families  of  curves  are  the  fit  obtained  with  B  «  5.17  K- 
cm^/J  with  the  Kirchhoff  transformation  (solid  lines)  and  B  »  4.05  K-cm^/J 
without  the  Kirchhoff  transformation  (dashed  lines)  using  a  » .0064  cm^/s,  the 
same  as  Fig.  2  for  the  following  full  widths  at  1/e  of  the  peak  intensity  of  the 
laser  beam  2(u  (mm):  2.03, 1.83, 1.65,  1.45,  1.24  (curves  from  left  to  right, 
respectively). 

(b)  TL  response  curves  from  Fig.  13  Ref.  5  recomputed  with  a  4. 18  W  Gaussian 
laser  beam  by  heating  TL  coating  through  the  glass  substrate.  The  dashed  curves 
are  obtained  with  B  «  3.35  K-cm^/J  and  the  solid  lines  are  curves  obtained  with 
B  ■  4.3  K-cm2/J,  using  the  Kirchhoff  transformation  and  a  ■  .0064  cm^/s  for 
the  following  measured  beam  widths  2q)  (mm);  2.13, 1.994, 1.65  (curves  from 
left  to  right,  respectively). 
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ABSTRACT 


The  theory  of  heating  a  semi-infinite  two-layer  system  with  a  localized  CO2  laser 
beam  of  uniform  circular  intensity  profile  is  applied  to  configurations  consisting  of  thin 
thermoluminescent  LiFrTiJvIg  phosphor  layers  on  borosilicate  glass  substrates.  We  study 
thermoluminescence  because  it  serves  as  a  convenient  monitor  for  the  temperature 
distribution  and  because  of  its  importance  in  developing  the  laser  heating  technique  for 
solid-state  dosimetry.  Experimental  and  calculated  results  are  compared  in  an  attempt  to 
completely  characterize  the  thermoluminescence  response  curves.  The  two-layer  system  is 
heated  in  two  different  modes:  the  laser  beam  impinges  onto  (a)  the  phosphor  layer,  and 
(b)  the  glass  substrate.  We  investigate  in  detail  changes  in  the  thermoluminescence 
response  due  to  different  laser  powers  and  variations  in  the  optical  and  thermal  properties 
of  the  samples. 


1.  INTRODUCTION 


bi  a  previous  publication  [1],  henceforth  lefened  to  as  I,  computational  methods 
were  presented  for  the  characterization  of  the  thermoluminescence  (TL)  response  obtained 
firom  thin  semi-infinite  thermoluminescent  phosphor  layers  on  transparent  substrates  upon 
exposure  to  localized  Gaussian  laser  beams.  These  luminescence  responses  (the 
thermoluminescence  intensity  measured  without  spatial  resolution)  are  not  the  familiar  TL 
glow  cu’^es  [2],  but  rather  exhibit  a  complex  shape  with  limited  glow  peak  separation, 
which  is  the  result  of  the  extreme  non-uniform  laser  beam  profile  and  the  resulting  non- 
uniform  temperature  distribution  in  the  sample.  At  any  given  time  after  onset  of  laser 
exposure,  the  lateral  temperature  distributi<m  is  neariy  Gaussian  as  well,  modified 
somewhat  by  thermal  diffusion.  The  spatially  resolved  TL  emission  from  the  thin  layer 
resembles  such  a  Gaussian  distribution  up  to  times  nearly  IJS  times  longer  than  the 
appearance  of  the  peak  of  the  spatially  unresolved  TL  response  curve  [3].  This  curve 
typically  exhibits  a  fast  rise  and  a  shallow  peak,  but  never  decays  completely  because  of 
lateral  heat  diffusion  [1].  The  TL  signal  disappears  only  when  the  laser  is  turned  off. 

Laser  heating  of  thin-layer  thermoluminescent  dosimeter  configurations  with 
focused  Gaussian  beams  has  been  applied  successfully  to  TL  imaging  [4]  and  to  rapid 
measurement  of  dose  distributions  generated  by  exposure  to  beta  rays  behind  attenuating 
niters  of  various  thicknesses  [5]. 

Fen*  accurate  and  sensitive  personnel  and  environmental  thermoluminescence 
dosimetry,  heating  with  a  laser  beam  of  uniform  intensity  profile  is,  however,  highly 
desirable.  As  we  will  demonstrate  below,  localized  laser  exposure  of  a  semi-infinite 
phosphor  layer  on  a  substrate  produces  glow  curves  which  are,  except  for  a  slowly 
decreasing  tail  at  the  high  temperature  side  of  the  peak,  nearly  identical  to  the  ones  obtained 
with  conventional  uniform  heating. 

In  this  paper  we  describe  a  TL  laser  heating  apparatus  that  features  a  circular  beam 
of  uniform  intension  profile  generated  with  a  technique  reported  in  Refs.  6  and  7.  It  was 
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employed  to  study  the  thermoluminescence  resptmse  curves  from  thin-layer  LiF:*Ii>lg 
dosimeters  on  glass  substrates  under  a  variety  of  experimental  conditions.  The  motivation 
for  this  investigation  was,  however,  not  limited  to  the  establishment  of  the  operating 
conditions  for  efficient  and  economical  laser  heating  in  thennoluminescence  dosimetry  and 
imaging  (the  latter  with  greatly  improved  spatial  resolution  as  compared  to  tire  non-uniform 
Gaussian  beam).  We  also  felt  that  the  thermoluminescence  emission  from  a  tiiin  layer  on  a 
substrate  is  a  very  convenient  monitor  for  the  time-evolution  of  the  temperature  distribution 
established  in  the  two-layer  system  during  locally  heating  with  a  laser  beam  of  uniform 
intensity  profUe.  The  general  theory  of  the  heat  transfer  from  such  a  beam  to  a  semi-infinite 
two-layer  system  and  the  subsequent  temperature  diffusion  was  presented  previously  [3,8]. 
In  the  following  we  will  frrst  briefly  review  the  basic  aspects  of  this  theory,  describe  the 
experimental  facility  some  detail,  and  finally  present  experimental  and  computation 
results. 

IL  COMPUTATIONAL  APPROACH 

In  applying  the  theory  of  Ref.  8  to  a  special  system  that  consists  of  a  thin 
thermoluminescent  phosphor  layer  on  a  strcmgly  absorbing  substrate,  a  number  of 
approximations  can  be  made  which  greatly  reduce  the  computational  effort  required  for 
simulating  thermoluminescence  response  curves  measured  under  various  different 
experimental  conditions.  The  justification  for  these  assumptions  is  discussed  in  detail  in 
Section  IV.  Here  we  only  briefly  review  the  impenant  features  of  the  simplified  theory. 

Assume  a  system  that  consists  of  a  thermoluminescent  phosphor  layer  on  a  glass 
substrate  in  perfect  thermal  contact  Each  layer  is  characterized  by  its  absorption  coefficient 
for  the  laser  photons,  thermal  conductivity,  Iq,  and  thickness  ly  The  index  i  s  1,2 
denotes  the  respective  layer  and  the  total  thickness  is  L  >  /j  /2.  In  contrast  to  the  general 

theory,  we  assume  for  the  special  case  of  a  thin  phosphor  layer  on  a  thicker  glass  slide  that 
the  thermal  diffusivity,  04,  is  temperature  independent  and  identical  for  both  layers.  The 


laser  beam  has  a  fixed  diameter  and  a  uniform  circular  intensity  profile  which  is 
cylindrically  symmetric  around  the  beam  axis.  It  enters  the  two-layer  system  perpendicular 
to  the  x-y  plane  either  from  the  side  of  the  phosphor  layer  or  the  glass  substrate.  The 
abso^ptance  (absorbed  fraction  of  the  incident  laser  power)  is  of  course  different  in  both 
cases,  and  so  is  the  co-called  piefactor  B  in  the  expression  fu’  the  tenq^erature  rise  (Eq.  1 
below).  The  thermal  diffusion  equation  and  the  apprc^riate  boundary  and  initial  conditions 
for  this  problem  are  stated  in  Ref.  8. 

Accounting  for  first-order  reflection  losses  on  the  surface  of  incidence  and  at  the 
layer  interface,  the  spatio-temporal  distribution  of  the  temperature  rise  for  the  two  layers  is 
given  by  [8] 

ATi(r,2,t)  =  (BP*/ti)2)  /o  d'C  Q(t-t){l  +  [(1  -  Ri)/JtB] 

0O 

2  exp(-^^/l)[exp(-n^/^)cos(Tl^,/l)  +  (Tlm/Jij)sin(Ti„/i)]/ 

I 

tNmd  +  llm^/4i2)cos(Tln^lll/l)]  +  (1  -  R2) 
exp(-Pi/j)(cos(iini/2)  +  (Tln/lia)sin(Tlm/2)  -  exp(-H2^2]/ 

[Nm(I  +Tlm2/^22)cos(Tlni/2)]}Xi(z)exp[-aTlm^(t-t)]  (1) 


where  P*  is  the  total  laser  power  (not  the  powCT  density)  incident  on  the  front  surface,  R 
the  reflectivity  layer  i,  (o  the  radius  of  the  uniform  circular  beam.  We  further  have 
B  -  (1  -  Ri)[l  -  (1  -  R2)  exp(-Hi/i  -  H2/2)  -  ^2  e*P('lAi/l)l/JtNo.  No  »  (1/a)  2  Iq/i ,  and 

X  l«l 

N„  ■  (1/2jc)  2  kj/iLl  +  sin20i/29i]  cos^Oj  with  9}  -.‘Hnifi.  and  Q(t)  ■  P(o))W2ot, 
thilax).  P(©4’)  -  /oto'do)'  Io(co',r)  exp[-((o'2  +  is  called  the  circular  coverage 
function.  The  parameters  Tim  ■")  ^  following  transcendental 

equation  (see  Appendix  A) 

ki  tan  Tim/i  + 1^2  tan  Tim/2  -  0,  (2) 


which  arises  from  the  boundary  conditions. 
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An  additional  simplification  of  Eq.  (1)  is  possible.  It  stems  from  the  fact  that  the 
material  of  interest  in  our  studies,  LiF,  is  not  a  strong  abso'ber  for  the  10.6  iim  photons  of 
the  laser  beating  beam  (p. «  40  cm'^).  In  contrast  to  this,  the  abstxpdon  coefficient  of  the 
glass  substrate  is  very  high.  Consequently,  the  laser  deposits  the  energy  directly  into  the 
interface  between  the  layers  if  it  impinges  cxi  the  phosphor  side  of  the  two-layer  system 
(Case  1  below).  On  the  other  hand,  it  is  possible  to  direct  the  laser  beam  onto  the  glass 
side  of  the  system  where  it  is  completely  absorbed  at  the  surface  (Case  2).  Thus  the  active 
thermoluminescent  layer  can  be  located  between  z  »  0  and  z-li  (Case  1)  or  between 

z  s  I,  to  z  s  L.  The  appropriate  expressions  for  the  temperature  rise  in  these  cases  are 
readily  obtained  from  Eq.  (1): 

Case  1:  thermoluminescent  layer  on  top  of  the  glass  substrate 

(4i/i  «  1) 

ATi(rAt)  =»  (BP*/co2)J^dTQ(r,t-t){l +No  I  cos(TimZ) 

«p[-«nni^  (t-X)]/[Nn,(l  +  cos(Tlro/i)] }  (3) 

Case  2:  thermoluminescent  layer  at  the  bottom  of  the  glass  substrate 
(42^2  «  1) 

AT2(rAt)  =  (BP*/<a2)Jo  dt  Q(r,t-x){l  +  Nq  I  cos(Tin,(L-z)) 

exp(-<*nm^  (t-X)]/[Nm(l  -t-  cos(Ti„/i)  cos(Tim/2)]}  (4) 

The  actual  computer  calculations  of  the  thermoluminescence  response  are  based  on 
the  well-known  first-order  electron  kinetics  of  the  LiF:Ti>Ig  phosphor  (TLD- 100, 
Harshaw/Filtrol )  and  require  knowledge  of  the  heating  rates  in  each  volume  element 
27irdidz.  These  are: 
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Case  1;  (p-i/i  «  1) 

aTitoVat  »  (BP*/(02)dtQ(r.t){l+No  I  cos(Tin^) 

exiH-anm^!l'[Nm(I  cos(Tln,/i)]}  (5) 

Case  2:  (^2^2  1) 

9T2(rAt)/3t  =  (BP*/o)2)dxQ(r,t){l +  No  S  cos(Tim(L-z)) 

msi 

expl-anm^«/INm(l  +  cos(ilm/i)  cos(Tln,/2)]}  (6) 

The  luminescence  glow  peaks  of  the  LiF*^g,Ti  considered  are  commonly  called  peaks  3, 

4,  and  5.  They  are  well  described  by  fint-order  kinetics  with  attempt-scape*frequency  of 
4  X  10^3  s-1^  7,3  X  1015  s*l,  and  4  x  10^1  s*l  and  trap  depth  of  1.23  eV,  1.54  eV,  and 
2.17  eV,  respectively  [9]. 

The  sum  of  the  relative  peak  heights  [9]  of  the  LiF;Ti>fg  glow  ctuve  is  used  as  a 
measure  for  the  total  trapped  electron  population,  and  the  relative  population  of  each  of  the 
three  electron  trap  levels  is  taken  to  be  proportional  to  its  corresponding  peak  height  with 
the  correction  of  luminescence  efficiency  at  different  heating  rates  [10].  TL  peaks  that 
occur  at  temperatures  lower  than  those  of  peaks  3  to  5  are  eliminated  by  a  preannealing  for 
10  minutes  at  100°C. 

The  three  first-order  electron-kinetic  differential  equations  for  peaks  3  to  S  [9]  are 
solved  with  these  heating  rates  by  the  Runge-Kutta  method  for  each  volume  element 
iTCtdcli  (for  case  1)  or  2}Crdr/2  (for  case  2).  Because  the  phosphor  layer  is  very  thin, 

integration  over  the  phosphor  layer  thickness  is  omitted  (see  Fig.  1 1  below).  In  Section  IV 
this  assumption  will  be  justified  by  comparison  of  measured  and  calculated 


thermoluminescence  response  curves.  The  numerical  integration  is  sufficiently  precise  with 
step  sizes  Ar  >  1.25oV50,  and  At »  10*^  tp,  where  tp  is  the  total  time  of  laser  exposure, 
except  for  the  first  millisecond  during  which  At «  10^  tp. 

in.  EXPERIMENTAL  SETUP 

The  experimental  arrangement  is  schematically  shown  in  Fig.  1.  A  temperature- 
stabilized  radio  frequency-excited  cw  CO2  laser  (Directed  Energy  Inc.  model  L3S-PS) 

emits  a  beam  whose  intensity  profile  has  an  annular  C'halo”)  cross  section  of  1  cm 
diameter.  The  10.6  pm  photon  beam  is  linearly  polarized.  The  maximum  output  power  is 
40  Watts.  For  closed-loop  control  of  the  laser  beam,  the  RF  power  supply  is  modified  for 
pulse-width  modulation  at  34  kHz  with  a  variable  duty  cycle  that  determines  the  laser 
output  power. 

The  beam  is  guided  to  a  focusing  module  by  two  parallel  gold-coated  copper 
minors.  Itentersthelight-tightopticalchamberthroughanAR-coatedGe  window.  The 
focusing  module  deflects  the  beam  down  onto  the  entrance  of  the  so-called  beam  equalizer 
channel  [6],  which  produces  a  uniform  square  beam  of  3  x  3  mm?  area.  To  simplify  the 
theory  and  the  numerical  calculation,  the  square  profile  of  the  output  beam  is  converted  to  a 
circular  one  with  the  aid  of  a  2.9S  mm  dianteter  aperture  at  the  channel  exit  The  sample 
holder  is  placed  in  close  proximity  to  it  so  that,  at  the  top  of  the  sample,  the  beam  has 
diverged  only  to  about  3  mm  diameter. 

A  gold-plated  electro-mechanical  shutto’,  placed  in  front  of  the  channel,  controls 
the  exposure  time  of  the  sample  to  the  laser  beam.  When  it  is  closed,  the  laser  beam  is 
reflected  into  a  beam  dump. 

To  provide  a  signal  for  closed-loop  control  of  the  laser  power,  an  anti-reflection 
coated  ZnSe  window  reflects  a  small  fixed  portion  (3  to  5  percent)  of  the  laser  beam  to  a 
piezoelectric  detector  (Baines  Engineering  Inc.  model  350-5  PZT).  A  beam  chopping 
wheel  operated  at  33.5  cycles/s  by  a  dc  motor  is  placed  in  front  of  the  detector  in  order  to 
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prevent  saturation.  The  ou^ut  of  the  detector  is  fed  to  the  electronic  control  unit,  providing 
an  updated  laser  power  measurement  every  9.3  ms.  In  this  way,  the  laser  output  power  is 
kept  constant  within  2-3  percent 

A  TL  samfde  on  a  25  x  25  mm?  substrate  can  be  placed  in  a  sample  holder  at  a 
distance  of 40-50  pm  from  the  end  of  the  charmeL  A  specific  area  of  the  sample,  selected 
by  a  pair  of  coupled  translation  stages,  can  be  heated  by  the  laser  beam.  In  this  way  up  to 
35  spots  on  a  continuous  TL  sample  can  be  heated  separately.  Typical  exposure  times 
range  from  100  ms  to  10  seconds.  The  translation  stages  are  connected  to  two  motcx' 
driven  actuators  which  are  controlled  by  a  motor  drive  system  with  a  positioning  accuracy 
of  up  to  0.1  piTL 

The  emitted  light  from  a  TL  sample  is  transmitted  by  a  light  guide,  through  a 
shutter,  and  a  protective  glass  window  to  a  photomultiplier  tube  (EMI  ^pe  9924B).  The 
shutter  and  glass  plate  prevent  the  PMT  from  ^xddental  exposxires  to  visible  light  or  the 
unimpeded  laser  power,  respectively.  The  light  guide  charmel  is  12.5  mm  diameter  buffed 
stainless  steel  tube.  To  maximize  the  collection  efficiency  of  the  emitted  thermo- 
luminescence,  one  end  of  it  is  placed  in  close  proximity  to  the  heated  area  on  the  sample. 
The  current  signal  of  the  PMT  is  measured,  processed,  stored,  and  displayed  with  the  aid 
of  a  transient  digitizer. 

The  electronic  control  unit  of  the  apparatus  performs  such  functions  as  timing, 
controlling,  and  protecting  the  system. 

When  a  TLD  reader  heating  cycle  is  activated,  the  following  sequence  is  performed 
by  the  control  electronic: 

The  laser  is  turned  on  and  its  output  power  is  measured  and  adjusted  until  the  laser 
is  stabilized  at  the  desired  power  level.  This  process  is  normally  accomplished  in  1(X)  to 
3(X)  ms;  however,  it  may  require  a  longer  time  if  the  laser  happens  to  be  in  an  undesirable 
operating  conditioiL  During  this  time,  the  shutter  remains  closed  and  the  sample  is  not 
exposed  to  the  beam.  After  the  desired  laser  power  is  obtained,  the  shutter  is  open  for  a  set 
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desired  heating  time.  At  the  end  of  the  heating  cycle,  the  shutter  will  close  and  the  laser  is 
turned  off. 

IV.  RESULTS  AND  DISCUSSION 

As  in  I,  we  have  chosen  LiF:Ti>fg  (TLD-100  supplied  by  Harshaw/Filtrol)  as  the 
thermoluminescent  material  to  serve  as  a  monitor  of  the  temperature  evolution  when  the 
two-layer  system  is  exposed  to  die  laser  beam.  All  systems  we  studied  consisted  of  a 
35-40  (im  thick  fine-grain  phosphor  layer  in  a  high-temperatuie  binder  on  150  ^  thick 
25  X  25  mm^  borosilicate  glass  microscope  cover  slides.  After  exposure  to  ^^o  gamma 
rays  (for  the  purpose  of  filling  the  traps)  and  a  subsequent  10  minute  pre-anneal  treatment 
to  eliminate  low-temperature  glow  peaks  [1],  the  center  of  these  samples  was  heated  with  a 
stable  and  uniform  3  ±  0.07  mm  diameter  CO2  laser  beam  whose  power  could  be  varied 

between  13J5  Watts  and  2.44  Watts.  The  laser  exposure  was  done  in  two  different  modes, 
the  laser  beam  either  striking  the  front  (phosphor  layer)  or  the  back  (substrate,.  Naturally, 
the  so-called  prefactor  B  (see  Section  II),  which  is  essentially  a  measure  of  the  absorptance, 
is  different  for  these  two  modes  of  laser  exposure.  With  the  method  described  in  I  its  value 
was  determined  to  be  B  >  10.7  ±  25%  K-cm^/J  for  the  phosph(»'  layer  surface  and 
B  s  6.3  ±  20%  K-cm^/J  for  the  glass  substrate  surface  (note  that  B  is  defined  slightly 
differently  as  compared  to  the  definidon  used  in  I;  here  it  is  divided  by  the  laser  power  P). 
The  electron-kinetic  data  of  the  traps  that  we  exploited  in  our  experiments  are  given  in 
Section  U. 

In  an  initial  set  of.  tests  we  established  the  validity  of  a  number  of  assumptions  and 
approximations  that  can  be  made  to  reduce  the  computational  effort  required  for  the 
comparison  of  experiments  and  the  theory,  without  severely  impacting  the  results. 

The  first  of  these  approximations  concerns  the  method  of  the  Kirchhoff 
transformation  described  in  I.  It  permits  one  to  account  for  the  temperature  dependence  of 
the  thermal  conductivity  of  the  glass  substrate,  which  we  approximated  as  k2  «  0.00942 


(1  -h  0.00222  AT2)  W/Kcm,  based  on  data  given  by  the  glass  manufacturer  [1 1].  In  Fig.  2 
the  TL  response  curve  is  shown  that  was  measured  by  exposing  the  glass  layer  (Case  2: 
heating  from  the  back  of  the  sample)  to  a  9.75  W  laser  beam  together  with  calculadois 
based  on  Eqs.  (4)  and  (6)  and  in  conjunction  with  the  simplifying  assumptions  discussed 
below.  The  KirchhofT  transformation  yielded  the  dashed  line,  and  the  solid  line  was 
obtained  assuming  constant  k2  >  0.00942  W/Kcm.  In  view  of  the  rather  good  agreement 

between  themy  and  experiment  that  is  achieved  without  accounting  fw  the  temperautre 
dependence  of  k,  we  decided  to  disregard  the  Kirchhoff  transformation  for  all  cases 
described  in  this  paper.  We  have  checked  this  with  a  number  of  different  experimental 
situations  and  found  it  was  always  justified  in  applications  of  the  theory  of  Ref.  8  to  thin 
two-layer  thermoluminescence  dosimeter  configuraticms. 

A  second  simplification  concerns  the  thermal  diffusivity  of  the  two  layers. 

Because  the  phosphor  layer  is  very  thin  and  does  not  strongly  absorb  the  laser  photons,  but 
the  glass  substrate  is  a  very  strong  absorber,  the  laser  energy  is  deposited  in  the  interface, 
when  the  system  is  heated  fiom  the  front  (Case  1).  The  phosphor  layer  is  indirectly  heated 
by  diffusion,  which  is  rather  rapid  in  z-direction,  because  its  so-called  thermal  response 
time  or  diffusion  time  [12],  tg  »  is  estimated  to  be  about  1  ms.  This  is  much  shorter 

than  the  time  required  to  measurably  change  the  thermoluminescence  intensity  even  at  the 
highest  laser  power  densities  that  appear  practical  for  this  type  of  heating  in  thermo¬ 
luminescence  dosimetry  [7].  Assuming  the  thermal  diffusivity  of  the  phosphor  layer  is  not 
significantly  greater  than  that  of  the  glass  substrate,  heat  transport  by  lateral  thermal 
diffusion  occurs  mainly  in  the  thicker  of  the  two  layen  with  essentially  instantaneous 
diffusion  in  z-direction  through  the  thin  layer.  This  is  of  course  also  true  for  Case  2  (laser 
heating  from  the  back).  Therefore,  we  are  justified  to  assume  the  thermal  diffusivity  of 
both  layers  to  be  equal,  its  value  being  that  of  glass  [3].  Actually,  there  is  a  reason  to 
assume  that  the  thermal  diffusivity  of  the  phosphor  layer  is  rather  small,  perhaps  even 
smaller  than  that  of  glass.  This  layer  is  a  loosely  packed  matrix  of  small  grains  in  a  poor 
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heat  conductor  (Dow  Coming  805  silicone).  Since  the  thermal  diffusivity,  oe,  for  glass  is 
slightly  temperature  dependent  (varying  from  0.0052  cm^/s  at  (y>C  to  0.0057°C  at  100°C, 
but  it  is  unknown  for  higher  temperatures),  and  since  we  did  not  choose  to  account  for  this, 
we  had  the  problem  of  selecting  a  proper  "temperature-averaged"  value  that  provided  a 
good  approximation  for  all  calculations  presented  below.  We  established  that  value  as  a  « 
0.0064  cmVs  by  the  best  Ht  of  families  of  TL  response  curves  measured  for  front  and  back 
exposure  with  several  different  laser  powers  (Bg.  3  is  an  example).  This  high  value 
appears  justified  when  one  considers  the  fact  that  the  TL  glow  peaks  we  selected  as 
"temperature  indicators"  for  our  experiments  reach  their  maximum  emission  intensity  at 
about  250°C  [9].  Consequently,  the  major  pcxtitm  of  the  thermoluminescence  emission 
and  thus  the  range  of  "temperature  monitoring"  occun  well  above  100°C.  It  should  be 
pointed  out  that  the  decay  of  glow  curve  is  very  soisidve  to  the  value  of  a  as  shown  in 
Fig.  4. 

Yet  another  simplification  of  the  computatitmal  procedure  is  possible  because  the 
TL  phosphor  layer  is  so  thin  and  diffusion  across  it .  •  relatively  rapid.  Instead  of 
calculating  the  transient  TL  emission  from  infmitesimally  thin  layers  5z  and  integrating 
over  z  from  z  «  0  to  z  « it  suffices  to  calculate  the  TL  emission  either  from  the  surface 
layer  at  z  «  0  or  from  the  interface  layer  at  z  « /j.  The  total  instantaneous  emission  is  then 
obtained  simply  by  multiplying  with  the  thickness,  /|,  of  the  active  layer.  In  each  situation 
good  agreement  results  between  theory  and  measurement  This  is  demonstrated  in  Bg.  5, 
where  the  case  of  Bg.  2  is  recalculated  for  different  values  of  z.  Henceforth,  we 
arbitrarily  choose  to  take  the  computed  TL  emission  from  the  interface  at  zs/^  to 

represent  the  emission  from  the  entire  layer. 

We  repeat,  these  simplifications  were  introduced  solely  to  reduce  the  computational 
effort,  because  the  so  executed  calculations  yield  results  that  agree  with  measurements 
within  the  experimental  errors.  They  are  not  strictly  required  as  the  more  exact  calculations 
pose  no  problems  in  principle. 


Even  though  we  have  not  yet  measured  in  detail  spatially  resolved 
thermoluminescent  resptmses  from  a  TL  phosphor  layer  on  a  glass  slide  that  is  heated  with 
a  localized  uniform  laser  beam,  we  present  here  some  initial  computadonal  results  as  diey 
help  to  visualize  the  reladon  between  the  spatio-temporal  temperature  distribution  and  the 
corresponding  thermoluminescent  response.  Figure  6(a)  shows  the  temperature 
distribution  in  the  phosphor-glass  interface  at  t «  4  ms  after  onset  of  exposure  of  the 
phosphor  layer  to  a  uniform  3  mm  diameter  laser  beam  of  6.5  Watt  power.  At  t » t*  >  1 15 
ms  (here  t*  is  the  time  of  TL  peak  occurrence),  the  distribution  is  considerably  different 
due  to  heat  diffusion  (Fig.  6b).  Interestingly,  the  l^eral  distributions  of  the  TL  emission 
appears  to  retain  the  basically  cylindrical  shape  up  to  t*;  it  simply  broadens  with  heating 
dme.  This  is  of  course  a  result  of  the  initial  rise  of  a  glow  curve  which  is  proportional  to 
exp(-E/kT)  (E  is  the  thermal  activation  energy).  Immediately  after  t*,  the  TL  intensity 
decreases  in  the  central  and  hottest  part  of  the  exposed  area,  and  where  it  has  decreased 
essentially  to  zero  at  t  *  U  t*,  shown  in  Figs.  7(a)-7(c).  For  times  t  >  t*,  the  TL 
emission  originates  only  from  those  areas  of  the  phosphor  layer  that  are  not  directly 
exposed  to  the  laser  beam,  but  are  heated  by  lateral  diffusion. 

For  applications  of  laser  heating  in  thermoluminescence  dosimetry  it  is  of  interest 
to  investigate  the  temperature  rise  of  the  exposed  phosphor  layer.  Optimum  operating 
conditions  are  commonly  believed  to  include  heating  a  thermoluminescent  dosimeter  with  a 
constant  heating  rate;  Le,  a  temperature  rise  linear  in  time  [13].  Figure  8  shows  computed 
temperature  rise  curves  for  the  same  sample  and  experimental  conditions  discussed  above. 
In  agreement  with  calculations  by  Sparks  [12],  the  temperature  in  the  center  of  the  exposed 
spot  rises  ftrst  noidinearly  for  a  few  ms  and  is  linear  thereafter.  The  heating  rate  is 
proportional  to  the  laser  power.  The  numerical  calculations  indicate  further,  that  the  heating 
rate  is  very  nearly  the  same  for  an  area  centered  around  the  laser  beam  axis  up  to  a  radius 
0  J  CO  for  heating  times  well  past  the  peak  of  the  TL  glow  curve  (compare  with  Figs.  4-7). 
This  is  of  considerable  practical  importance  as  it  demonstrates  by  the  usefulness  of  laser 


heating  not  only  for  dosimetry  [3, 14, 15],  but  also  in  a  number  of  special  areas  which 
exploit  thermoluminescence;  e.g.,  dating  of  archaeological  artifacts  and  art  authentication 
[IQ.  Very  precise  Unear  and  reproducible  heating  is  possible  of  minute  fine-grain  samples 
diat  are  attached  with  a  suitable  binder  to  a  substrate  in  the  center  of  the  uniform  laser 
heating  beam.  We  have  successfully  heated  single  <100  |jm  diameter  grains  of  natural 
quartz  and  obtained  strong  TL  signals  without  previous  additional  radiation  exposures  [17]. 
The  reproducibility  and  accuracy  of  this  heating  technique  depends,  of  course,  on  the 
stability  and  reproducibiUty  of  the  laser  beam.  As  mentioned  in  Section  HI,  it  is  possible  to 
control  the  laser  power  to  within  less  than  3  percent 

One  of  the  goals  of  our  investigations  was  to  establish  the  vaUdity  of  the  general 
laser  heating  theory  presented  in  Ref.  8.  An  assessment  of  the  achievable  accuracy  of  the 
thermcluminescence  temperature  monitoring  method  was  therefore  essential.  Experimental 
errors  arise  from  the  uncertainty  in  laser  beam  parameters,  the  phosphor  layer  uniformity, 
and  the  reproducibiUty  with  which  the  two-layer  samples  can  be  produced.  With  a  stable 
laser  power,  only  the  beam  radius,  w,  remained  as  a  critical  parameter.  Calculations 
indicated  that  it  very  sensitively  influences  the  TL  glow  curves.  Changing  O)  by  only  ±1 
percent  yields  shifts  (Fig.  9)  that  otherwise  are  produced  only  by  variation  of  the  produced 
layer  thicknesses  (see  below).  For  this  reason  the  laser-heated  thermoluminescence 
apparatus  described  in  Section  HI  was  constructed  with  a  fixed  rather  than  a  variable  laser 
beam  radius,  and,  therefore  no  experiments  were  performed  to  date  with  a  heating  beam  of 
constant  power  but  variable  diameter. 

Variations  in  the  thickness,  /j,  of  the  phosphor  layer  cause  the  TL  response  curves 
to  shift  and  change  in  intensity.  This  is  demonstrated  in  Figs.  10  and  11,  which  show  a 
number  of  curves  that  were  measured  (Fig.  10)  and  calculated  (Fig.  1 1)  with  different 
values  of  l\.  Figure  10  indicates  the  rather  lintited  reproducibility  of  the  spreading  method 

used  to  fabricate  the  samples.  Despite  considerable  effort,  which  involved  the  construction 
of  a  special  layer  spreading  machine  and  careful  control  of  the  phosphor-binder  mixture, 


this  performance  could  not  be  improved  beyond  the  one  reported  here,  and  more  work  on 
layer  fabrication  techniques  is  needed.  However,  the  layer  uniformity  achieved  was 
sufficient  for  a  study  of  the  TL  response  of  a  number  of  nominally  identical  samples  as  a 
functicm  of  laser  power.  Typical  measured  and  computed  glow  curves  resulting  from  this 
study  are  shown  in  Fig.  12  (Case  1:  the  laser  heating  beam  impinges  onto  the  TL  layer). 
For  a  best  fit,  the  parameter  B  had  to  be  adjusted  slightly.  Note  that  the  variation  in  B  is 
within  the  accuracy  of  the  experimental  method  we  used  to  determine  it  [1].  It  was 
necessary  to  also  account  for  the  fact  that  the  overall  luminescent  yield  decreases  with 
increasing  heating  rate  for  otherwise  idential  expoimental  conditions.  Detailed 
measurement  of  the  dependence  of  the  yield  on  the  heating  rate  will  be  presented  in  a 
forthcoming  publication  [10]. 

V.  CONCLUDING  REMARKS 

We  have  compared  experimental  and  computed  thermoluminescence  response 
curves  obtained  from  two-layer  configurations  by  localized  heating  with  a  3  mm  diameter 
CO2  laser  beam  of  uniform  circular  intensity  profile  of  power  up  to  13  Watts.  The  results 

confirm  the  theory  of  Ref.  8  and  show  that  thin-layer  thermoluminescence  dosimeters  can 
be  heated  with  heating  rates  up  to  SSOCPOs. 

It  is  of  particular  importance  that  the  thermoluminescence  glow  curves  measured 
under  these  experimental  conditions  have  the  familiar  shape  obtained  with  conventional 
contact  heating. 
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APPENDIX  A 

In  this  work  we  have  a  poorly  absorbing  phosphor  coating  /j  «  0.2L  on  a  highly 
absorbing  glass  slide,  I2  •  0.8L  where  the  total  thickness  L  *  0.01875  cm.  In  this  case  m 
distinct  toots  1101  can  be  obtained  from  the  transcendental  equation  k^  tan  0 
kj  tan  40-0  with  0-O.2ti,„L.  With  F  -  5jc/L,  h  -  3+2Ac,  k  -  kj/kj,  D+ -  {h  +  [h^ - 
(1  +  and  D*  -  {h  -  [h  +  (1  +  4/k)]^^}^,  then  the  m  roots  greater  than 

Tio-O  are  given  by 

Tim  •  1’'“^  tan-i(D-)  +  (m  - 1)]  •  F  for  m  -  1,  6,  1 1,... 

[tr^  tan**(D+)  +  (m  - 1)]  •  F  for  m  -  2,  7,  12,... 

[1  -  jr^  tan-i(D*)  +  (m  -  1)]  •  F  for  m  -  3,  8, 13,... 

[1  -  tt'  tan*i(D+)  +  (m  •  1)]  •  F  for  m  »  4,  9,  14,... 

m-  F  for  m  »  5, 10,  15,... 

The  root  -  0  yields  the  first  tenn,  unity,  in  the  brackets  of  Eq.  (1). 
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FIGURE  CAPTIONS 

Fig.  1 .  Schematic  experimental  arrangement  of  the  laser  theimoluminescence  apparatus. 

Fig.  2.  Thermoluminescence  response  curves  measured  and  calculated  with  a  uniform 
laser  beam  of  3  mm  diameter  and  9.75  Watts  power  heating  a  37  J  nm  thick  TL 
layer  through  a  0.15  mm  thick  b<»osilicate  glass  substrate  (the  laser  beam 
impinges  onto  the  glass  substrate).  The  solid  line  is  computed  with  a  constant 
thermal  conductivity  of  the  glass  k  »  0.(X)942  W/Kcm  and  the  dashed  line  with 
the  method  of  the  Kirchhoff  transformation,  which  takes  into  account  the 
temperature  dependence  of  k[ll.  With  the  thermal  diffusivity  a  =  0.0064  cm^/s 
(see  Fig.  3),  the  parameter  B  was  determined  to  be  6.4  Kcm^/J  by  the  best  fit  to 
the  measured  curve.  This  value  agrees  well  with  the  experimentally  determined 
one:  B  ■  5.7  ±  20%  Kcm^/L  (see  Ref.  1).  Prior  to  the  measurement  of  the  TL 
response  curve,  the  sample  was  exposed  to  gamma  rays  from  a  ^o  source  and 
pre-annealed  for  10  minutes  at  100°C. 
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Fig.  3.  Thermoluminescence  response  curves  measured  by  headng  with  laser  beams  of 
fixed  3  mm  diameter  but  different  powers  impinging  onto  the  glass  side  of  the 
two-layer  system.  The  sample  is  a  37.5  pm  thick  liFtH^g  on  a  0. 15  ram  thick 
boiosilicate  glass  substrate.  The  family  of  curves,  calculated  from  Eqs.  (4)  and 
(6),  yields  the  value  a  «  0.QQ64  cm^/s  for  the  thermal  diffusivity  via  best  fit  with 
the  measurements.  The  calculations  were  performed  with  k^ »  0.00942  W/Kcm 
that  is  characteristic  for  the  system  and  the  following  laser  powers  and  B 
parameters: 


13 

Watts 

B»  6.3 

Kcm^/J 

9.75 

Watts 

6.4 

Kcm^/J 

6.5 

Walts 

6.0 

Kcm^/J 

4,875 

Walts 

6.3 

Kcm^/J 

Peaks  I  and  2  of  LiF:Ti>lg  (TLD-100,  Harshaw/Filtrol)  were  removed  by  10 
minutes  preaimealing  at  lOOX  [9]. 

Fig.  4.  Computed  thermoluminescence  response  curves  of  the  sample  from  Fig.  3  for 

various  thermal  diffusivities,  ol  The  best  fit  with  the  measured  curve  is  obtained 
with  a  »  0.0064  cm2/s  and  the  prcfactor  B  *  6.0  Kcm^/J  (dotted  curve).  The 
effect  of  different  thermal  diffusivities  is  demonstrated  by  the  solid  (a  =  0.005 
cm2/s)  and  the  dashed  (a  *  0.0078  cra^/s)  lines.  Laser  power:  P*=  6.5  Watts. 


Fig.  5.  Thennoluininescence  response  curves  obtained  for  the  sample  and  experimental 
conditions  as  in  Fig.  4,  computed  for  different  depths,  z,  of  the  active  layer.  The 
dashed  line  represents  the  normalized  tlwrmoluminescence  emissi  jn  from  a  thin 
sheet  at  z  3  L  (the  surface  of  the  TL  layer),  and  the  solid  line  that  at  z  » (the 

interface  of  the  TL  layer  and  the  glass  substrate).  We  have  adjusted  the  parameter 
B  so  that  the  occurrence  of  the  peak  emission  is  the  same  as  that  of  the  measured 
curve:  B  =  6  Kcm^/J  (dashed  line)  and  B  »  5.2  Kcm^/K  (solid  line).  The 
calculations  demonstrate  that  for  thin  TL  layers  the  z  dependence  of  the  TL 
response  is  small. 

Fig.  6.  Computer  simulations  for  the  normalized  spatial  temperature  distributions  T(r,z,t) 
calculated  for  exposure  of  a  phosphor  layer  on  a  glass  substrate  to  a  6  J  W-laser 
heating  beam  with  a  uniform  circular  intensity  profile,  2w  »  3  mm.  The 
calculations  are  performed  with  B  »  7  Kcm^/J,  a  »  0.(X)64  cm^/s  and  a  total 
thickness  L  =  0.1875  mm  for  various  times  after  onset  of  laser  exposure:  a)  t  =  4 
ms  and  b)  t « t*  *  1 15  ms  (the  peak  intensity  of  the  thermoluminescence  response 
occurs  at  t). 

Fig.  7.  The  normalized  spatial  distributions  of  the  thermoluminescence  response 

calculated  for  the  sample  from  Fig.  6  at  various  times  with  respect  to  the  time,  T’*', 
of  the  TL  peak  intensity  occurrence,  t  which  are:  a)  t  =  1/2 1^,  b)  t  =  t*,  and 
c)t»1.5t*. 

Fig.  8.  Computed  time  evolution  of  the  temperature  of  the  sample  and  experimental 

conditions  as  in  Fig.  6  at  various  distances,  r,  from  the  laser  beam  axis:  r  »  0, 
«i/4,  w/2,  w,  1. 125  w  (from  top  to  bottom),  (w  *  0. 15  cm  is  the  radius  of  the 
laser  beam.) 
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Fig.  9.  Calculated  thermoluminescence  response  curves  of  the  sample  from  Fig.  3  for 
small  variation  in  the  radius  of  the  laser  beam.  The  computations  are  performed 
with  P*  «  6  J  Watts  and  the  following  beam  radii  w  (from  left  to  right): 

0.148S,  0.15,  and  0.1515  cm. 

Fig.  10.  A  set  of  thermoluminescence  response  curves  measured  with  a  37. 5  (im  thick  TL 
layer  on  0.15  mm  thick  borosilicate  glass  substrate  for  fixed  laser  beam  radius  w 
a  0.15  cm  and  various  laser  powers  P*.  The  laser  beam  impinges  on  the 
phosphor  side  of  the  sample.  To  check  the  reproducibiliQ'  of  the  measurements, 
several  thermoluminescence  responses  fiom  different  sites  on  the  sample  were 
measured  for  otherwise  identical  conditions  with  the  following  laser  powers  (in 
Watts):  1)  13.  2)  9.75,  3)  6.5. 4)  4.875, 5)  3.25,  and  6)  2.44. 

Fig.  11.  Computed  normalized  thermoluminescence  response  curves  for  various  total 

thicknesses,  L,  but  fixed  ratio  IxU'i  =  1/4.  Otherwise  the  sample  is  identical  to  that 

of  Gi.  3.  The  calculations  were  performed  with  a  laser  power  of  6.5  Watts  and 
the  following  thicknesses  (fi-om  left  to  right):  0.094, 0.178, 0.1875, 0.197,  and 
0.375  mm. 

Fig.  12.  Measured  and  computed  thermoluminescence  response  curves  of  the  sample  from 
Fig.  10,  and  the  same  experimental  conditions.  The  dashed  curves  are  computed 
fiom  Eqs.  (3)  and  (5)  and  the  first-order  electron  kinetics  parameters  determined 
by  McKeever  [9]  for  peaks  3-5  of  LiF:Ti>Ig  and  by  taking  into  account  the 
measured  decrease  in  luminous  efficiencies  with  increasing  heating  rates.  Again, 
peak  2  was  removed  by  a  10  minute  preaimealing  at  lOOX.  With  the  thermal 
diffusivity  a  >  0.0064  cm^/s  and  a  laser  beam  radius  of  w  s  0.15  cm,  the 
parameter  B  was  determined  for  each  curve  by  the  best  fit  The  results  are 
calculated  with  the  following  luminescence  efficiencies,  ii(in  %  of  the  value 
obtained  for  the  smallest  heating  rate):  81.8,  81.9,  85.5,  89.6,  95.7,  and  100 
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DEVELOPMENT  OF  OPTICALLY  STIMULABLE  STORAGE  PHOSPHORS 
FOR  IMAGING  OF  NEUTRON  DOSE  DISTRIBUTIONS 

I.  INTRODUCTION 

In  studies  that  commenced  in  1981  we  have,  in  collaboration  with  Prof.  J. 
Gasiot,  University  of  Montpellier,  France,  established  that  certain  optically 
stimulable  phosphors  can  retain  a  fraction  of  the  energy  absorbed  from 
ionizing  radiation  fields  such  as  UV,  x-rays,  protons,  etc.  and  that  this 
energy  can  be  released  in  form  of  visible  luminescence  on  a  time  scale  of 
about  100  ns  when  exposed  to  an  infrared  laser  beam.  Furthermore,  the  infrared 
stimulation  spectrum  of  some  of  these  phosphors  turned  out  to  overlap  the 
emission  wavelength  of  conventional  Nd:YAG  lasers.  This  finding  pointed 
Lmmediately  to  applications  in  real-time  imaging  of  dose  distribution  [1],  the 
most  intriguing  of  which  appeared  to  be  neutron  imaging  for  the  purpose  of 
non-destructive  testing  of  large  objects  containing  elements  of  low  order 
number,  e.g.  polymers.  The  mode  of  operation  one  can  envisage  is  the  use  of  a 
so-called  proton  radiator  in  form  of  an  approximately  1  mm  thick  polyethylene 
layer  in  front  of  the  phosphor  screen.  The  neutrons,  generated  in  a 
conventional  neutron  generator  rather  than  a  reactor  or  a  neutron  source  such 
as  Pu-Be  or  252  Cf  (to  keep  the  gamma  component  of  the  radiation  as  low  as 
possible),  produce  knock-on  protons  in  the  proton  radiator  which  in  turn 
interact  with  the  grains  of  the  phosphor.  In  this  process  the  electron  traps 
are  filled  proportionally  to  the  proton  density.  Thus,  a  latent  image  is 
stored  which  is  released  upon  exposure  to  the  stimulating  IR  laser  beam.  This 
beam  will  scan  the  phosphor  and  in  this  way  produce  a  luminescence  emission 
that  is  detected  by  a  photomultiplier  tube,  digitized,  and  processed  by  a 
commercially  available  imaging  processor.  In  contrast  to  other  luminescence 
converter  phosphors,  the  storage  capability  of  the  phosphors  we  investigate 
facilitates  a  mode  of  operation  that  might  be  called  real-time  or  nearly 
real-time  imaging  with  intermediate  information  up-dating:  while  a  pixel  or  a 
part  of  phosphor  screen  is  "read”  by  the  scanning  IR  laser  beam,  the 


radiation  can  already  imprint”  new  information  on  the  section  of  the  screen 
that  is  not  exposed  to  the  stimulating  laser  reading  beam.  The  scanning  laser 
beam  is  an  ideal  way  to  digitize  the  image  on  a  pixel-  by-pixel  basis.  For 
real-time  scanning-type  digitization  to  be  possible,  no  luminescence 
overlapping  the  IR  stimulated  luminescence  band  from  a  pixel  exposed  to  the 
IR  laser  beam  can  be  allowed  to  be  emitted  from  the  entire  screen.  This  is 
particulartly  true  for  the  luminescenece  emitted  during  exposure  to  the 
ionizing  radiation.  In  other  words,  the  luminescence  emitted  during  and  after 
exposure  of  the  screen  to  the  ionizing  radiation  must  occur  in  an  emission 
band  that  is  different  from  the  IR  stimulated  emission  so  that  it  can  be 
completely  blocked  with  filters  from  reaching  the  detector. 

During  examination  of  a  large  number  of  IR  stimulable  phosphors,  mostly 
prepared  and  provided  Prof.  Gasiot  for  our  use  prior  to  this  contract,  we 
observed  that  certain  samarium-cerium  co-doped  sulphides  emit  orange 
luminescence  during  UV  excitation,  but  switch  to  a  distincly  green  emission 
when  exposed  afterwards  to  IR.  This  more  or  less  accidental  observation  formed 
the  basis  for  the  present  study  whose  goal  is  to  develop  the  phosphors  for 
later  use  in  nearly  real-time  neutron  imaging. 

Of  course,  the  phosphor  must  posses  certain  properties  to  be  a  successful 
candidate  for  the  future  development  of  a  complete  neutron  imaging  system 
comprising  a  long  "wish  list”.  It  is  the  purpose  of  the  work  performed  under 
this  contract  to  identify  such  suitable  phosphors  and  to  establish  their  basic 
physical  characteristics  vis-a-vis  the  goal  of  near  real-time  neutron  imaging. 

The  ideal  phosphor  must: 

1.  emit  a  different  luminescence  spectrum  during  exposure  to  the 
radiation  whose  distribution  is  to  be  imaged  as  compared  to  the  spectrum 
emitted  during  image  retrieval  by  the  scanning  IR  laser  beam, 

2.  be  efficiently  stimulable  by  a  conventional  commercially  available 
laser  whose  wavelength  must  be  such  as  to  not  cause  a  signal  in  the  detector 


1% 


used  for  signal  generation  with  the  IR  stimulated  luminescence,  i.e.  the  IR 
stimulation  spectrum  must  be  different  from  the  IR-stimulated  luminescence 
emission  spectrum, 

3.  be  capable  of  retaining  an  acceptable  fraction  of  the  energy  from 
the  ionizing  radiation  (most  importantly  neutron-generated  knock-on  protons), 
and  efficiently  convert  the  stored  energy  to  luminecence  upon  stimulation  by 
the  IR  laser  beam, 

4,  be  capable  of  a  large  dynamic  range,  that  is  the  dose 
proportionality  must  have  a  useful  range  of  at  least  four  orders  of  magnitude. 

While  it  is  well  known  that  the  conversion  efficiency  of  energy  absorbed  from 
the  radiation  field  to  visible  luminescence  can  approach  several  percent  and 
the  IR  stimulation  spectra  and  emission  spectra  of  a  number  of  phosphors  have 
the  desired  properties  listed  above  under  2.,  producing  a  phosphor  that  emits 
a  different  luminescence  spectrum  during  exposure  to  the  ionizing  radiations 
as  compared  to  the  one  generated  during  IR  stimulation  presented  a 
considerable  challenge.  As  this  property  is  the  most  important  one,  we  have 
concentrated  on  it  during  the  contract  period  from  7-15-85  through  3-15-86. 

The  phosphors  we  have  selected  for  this  study  are  MgS;Ce(0. 1%) .  Sm(0.05%)  and 
CaS;Ce(0. 1%) ,Sm(0.05S) ,  They  represnt  the  best  phosphors  prepared  to-date. 

In  the  following  we  will  first  decribe  the  method  of  their  preparation,  which 
has  emerged  in  extensive  trials  with  differnt  preparation  parameters. 
Thereafter,  we  report  on  the  spectral,  decay  and  storage  properties  of  the 
luminescence. 
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beam  is  an  ideal  way  to  digitize  the  image  on  a  pixel-  by-pixel  basis.  For 
real-time  scanning-type  digitization  to  be  possible,  no  luminescence 
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During  examination  of  a  large  number  of  IR  stinulable  phosphors,  mostly 
prepared  and  provided  Prof.  Gasiot  for  our  use  prior  to  this  contract,  we 
observed  that  certain  samarium-cerium  co-<loped  sulphides  emit  orange 
luminescence  during  UV  excitation,  but  switch  to  a  distincly  green  emission 
when  exposed  afterwards  to  IR.  This  more  or  less  accidental  observation  formed 
the  basis  for  the  present  study  whose  goal  is  to  develop  the  phosphors  for 
later  use  in  nearly  real-time  neutron  imaging. 

Of  course,  the  phosphor  must  posses  certain  properties  to  be  a  successful 
candidate  for  the  future  development  of  a  complete  neutron  imaging  system 
comprising  a  long  "wish  list".  It  is  the  purpose  of  the  work  performed  under 
this  contract  to  identify  such  suitable  phosphors  and  to  establish  their  basic 
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1.  emit  a  different  luminescence  spectrum  during  exposure  to  the 
radiation  whose  distribution  is  to  be  imaged  as  compared  to  the  spectrum 
emitted  during  image  retrieval  by  the  scanning  IR  laser  beam, 

2.  be  efficiently  stimulable  by  a  conventional  commercially  available 
laser  whose  wavelength  must  be  such  as  to  not  cause  a  signal  in  the  detector 


2.  CaF  sub2  and  KBr  are  inefficient  fluxes. 

3.  NaCL  and  NH  sub  4  Cl  are  more  efficient  than  MgCL  sub  2  and  BaCl 
sub  2,  but  all  work  fine. 

4.  Increasing  the  flux  concentration  improves  IR  stimulability  while  simulta¬ 
neously  decreasing  the  Sm  luminous  efficiency. 

5.  The  most  effective  flux  concentrations  turned  out  to  be  in  the  range  from 
0.5%  to  4.0%. 

III.  LUMINESCENCE  PROPERTIES  OF  CE/SM  CO-DOPED  MAGNESIUM  SULPHIDE  AND  CALCIIT-I 
SULPHIDE 

The  properties  investigated  to-date  include  the  IR  stimulation  spectra,  the 
emission  spectra  during  UV  excitation,  during  simultaneous  UV  excitation  and 
IR  stimulation,  and  during  IR  stimulation  at  various  storage  times  after  UV 
excitation.  All  these  experiments  were  carried  out  with  the  aid  of  a  Princeton 
Applied  Research  optical  multichannel  analyser  (type:  OMA  III). 

Initial  experiments  have  also  been  performed  on  the  decay  times  of  the 
luminescence  emission,  but  a  detailed  investigation  of  the  IR  stimulation 
properties  (response  time,  decay  time  as  a  function  of  power  density,  and  the 
experimental  conditions  required  to  completely  empty  the  traps  and,  thus, 
quench  the  stored  image  information  by  IR  stimulation)  has  not  yet  been 
executed.  It  requires  the  use  of  a  NdrYAG  laser  system,  which  will  be 
available  for  these  experiments  after  May  15,  1985. 

1.  Luminescence  Emission  Spectrum  during  UV  Excitation 

Exposure  of  the  phosphor  layers  to  the  253.7  nm  mercury  line  (mercury  penlight 
plus  narrow  band  interference  filter)  yields  the  emission  spectra  of  Figs.  1 
and  2.  Both  exhibit  the  characteristic  lines  of  Sm  and  a  broad  unresolved  band 
due  to  Ce  [2] . 
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II.  PHOSPHOR  PREPARxMION 

Reproducible  results  in  preparing  IR  stimulable  phosphors  can  only  be  achieved 
by  careful  process  control.  Our  preparation  facility  consistes  of  a  Lindberg 
three  zone  tube  furnace  and  specially  designed  gas  handling  equipment. 

Phosphor  preparation  includes  the  following  steps: 

1.  Reduction  of  MgS04  or  CaS04  to  MgS  or  CaS  in  an  atmosphere  of  dry  nitrogen 
and  carbon  disulphide.  Reagent  grade  nitrogen  is  dried  by  bubbling  through 
fuming  sulfuric  acid.  Carbon  disulfide  is  admixed  by  bubbling  this  dried  gas 
through  liquid  CS  sub  2.  This  gas  mixture  enters  the  aluminum  oxide  tube  in 
the  furnace,  in  the  center  of  which  the  powders  to  be  reduced  are  placed  in 
quartz  boats  or  aluminum  oxide  containers.  Quartz  permits  faster  cooling 
times,  a  feature  necessary  for  faster  preparation  of  MgS  phosphors.  The 
reduction  process  requires  heating  the  material  in  the  CS  sub  2  /N  sub  2 
atmosphere  at  850*C  for  one  hour. 

2.  The  sulphides  produced  in  step  1.  are  intimately  mixed  with  the  dopants  (in 
form  of  samarium  oxide  and  cerium  oxide),  the  flux  agent  (see  figures  for 
details),  and  additional  sulfur. 


3.  This  mixture  is  fired  in  dry  nitrogen  gas  at  1050®C  for  MgS  and  at  1220®C 
for  CaS.  After  cooling  to  SSO^C,  the  CS  sub  2  gas  is  added  to  the  nitrogen. 
The  total  firing  cycle  requires  two  hours  for  proper  diffusion  of  the  dopants 
through  the  powder  grains. 

To  produce  these  best-effort  phosphors,  carefully  searching  for  the  optimum 
heating  cycles  was  required.  In  addition,  we  had  to  experiment  with  various 
concentrations  of  the  dopants  as  well  as  different  fluxing  agents  of  various 
concentrations.  Fluxes  tested  included  CaF  sub  2,  KBr,  ;IH  sub  4C1,  NaCl,  MgCl 
sub  2,  and  3aCl  sub  2.  These  flux  tests  can  be  summarized  as  follows; 

1.  Sm  requires  no  flux  for  efficient  diffusion. 
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In  view  of  the  anticipated  application  of  both  phosphors  in  radiation  imaging 
and  the  requirement  that  IR  stimulation  must  be  possible  with  a  cost-efficient 
scanning-beam  laser  system  such  as  a  Nd:YAG  laser,  Che  range  of  IR  wavelengths 
chat  can  stimulate  the  luminescence  emission  by  releasing  the  electrons  from 
the  storage  traps  is  of  particular  Interest.  We  have  measured  these  spectra 
with  an  IR  monochromator  (gracing  blazed  for  1.3  urn).  The  results,  still 
preliminary  in  detail(e.g.  the  apparent  structure  is  unconfirmed)  are 
presented  in  Figs.  7  and  8.  As  hoped,  the  stimulation  spectra  for  both 
phosphors  peak  at  or  near  the  1.06  um  wavelength  of  the  )ld;YAG  laser  photons. 


5.  Decay  and  Storage  Properties 


Since  admixtures  of  both  spectral  components  in  the  UV  excitation  and  IR 
stimulation  spectra  are  always  present  and  complete  spectral  discrimination  in 
Che  excitation  and  stimulation  mode  of  the  imaging  process  does  not  appear 
feasible  with  these  phosphors,  we  have  investigated  the  possibility  of  near 
real-time  imaging.  Here,  the  excitation  with  knock-on  protons  or  other 
ionizing  radiation  does  not  occur  during  IR  stimulation  v/ich  the  scanning 
laser  beam.  Instead,  the  excitation  source  (i.e.  the  neutron  generator  beam) 
is  interrupted  during  the  scanning  process.  Since  we  fully  expect  the  response 
time  to  the  scanning  beam  to  be  very  fast,  rapid  "reading"  of  the  stored  image 
with  a  sufficiently  intense  Nd:YAG  laser  beam  [1]  is  possible.  Thus,  the  image 
can  be  updated  in  as  little  time  as  it  Cakes  to  store  a  new  image  (which 
depends  of  course  on  the  neutron  fluence)  plus  Che  times  required  for  the 
complete  decay  of  Che  undesired  spectral  component  present  during  UV 
excitation  (in  our  case  the  Ce  band)  and  for  reading  with  the  scanning  laser 
beam.  As  was  shown  previously,  the  reading  time  can  be  much  less  than  a  second 
[ 1 ] .  We  expect  the  imaging  system  to  provide  an  updated  image  in  less  then  a 
few  seconds. 

Therefore,  knowledge  of  the  decay  time  of  the  undesired  spectral  component, 
the  Ge  band  present  during  UV  stimulation,  is  of  paramount  importance.  We  have 
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2.  Luminescence  Emission  Spectra  during  Simultaneous  UV  Excitation  and  IR 


Stimulation 


An  indication  of  the  effectiveness  of  IR  stimulation  may  be  had  by 
simultaneously  exposing  the  phosphor  to  UV  and  IR.  The  IR  illumination  was 
obtained  from  a  Kodak  slide  projector  (300  Watt  quartx-halogen  Lamp,  heat 
filter  replaced  by  2  mm  thick  polished  but  uncoated  silicon  window  to  block 
all  radiation  below  900  nm  wavelength).  Figs.  3  and  4  demonstrate  that  IR 
exposure  in  addition  to  253.7  nm  UV  excitation  significantly  enhances  the  Ce 
emission  band,  but  does  not  affect  the  Sm  lines. 


3.  IR  Stimulated  Luminescence  Emission  Spectra 


After  filling  the  electron  storage  traps  by  UV  excitation,  exposure  to  IR 
causes  luminescence  emission  whose  spectra  are  shown  in  Figs.  3  and  6.  They 
clearly  demonstrate  that  the  orange  component  that  is  so  dominant  during  UV 
excitation  is  now  very  weak  compared  to  the  green  Ce  emission  band.  This  is 
precisely  the  hoped  for  property  of  the  phosphors  we  have  developed  for 
eventual  later  application  in  radiation  imaging.  In  comparing  the  UV  excited 
luminescence  with  the  IR  stimulated  one,  it  also  becomes  clear  that  the 
anticipated  complete  separation  of  these  spectra  is  not  present.  Both  types 
always  contain  components  of  both  the  Sm  and  the  Ce  emission  spectra.  Complete 
separation  would  have  meant  that  real-time  imaging  is  feasible.  In  light  of 
the  admixture  of  the  Sm  lines  in  the  IR  stimulated  spectra  and  of  the  Ce  band 
in  the  UV  excited  ones,  only  "nearly  real-time"  imaging  appears  possible.  We 
will  discuss  this  point  in  more  detail  below. 


4.  IR  Stimulation  Sectra 
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FIGURE  CAPTIONS 


Fig.  1 

Luminescence  emission  spectrum  for  MgS:Ce(0.1%) ,Sm(0.05%)-t-2%  Cl  flux  [standard 
MgS  phosphor]  during  excitation  with  253.7  nm  UV.  Spectral  assignment  from 
Ref.l 

Fig.2 

Luminescence  emission  spectrum  for  CaS:Ce(0.1%)  ,Sm  0.05%)-t>2%  Cl  flux  [standard 
CaS  phosphor]  during  excitation  with  253.7  nm  UV.  Spectral  assignment  as  in 
Fig.l 

Fig. 3 

Luminescence  emission  from  the  standard  MgS  phosphor  obtained  by  simultaneous 
excitation  with  253.7  nm  UV  and  stimulation  with  broadband  IR  (quartz  halogen 
lamp  plus  5  mm  thick  Si  filter) 

Fig.  4 

Same  as  Fig.  3  for  standard  CaS  phosphor 
Fig.  5 

Luminescence  emission  spectrum  from  the  standard  MgS  phosphor  after  excitation 
with  253.7  nm  UV  and  during  subsequent  stimulation  with  broadband  IR.  Mote  the 
difference  compared  to  Figs.  1  and  3 

Fig. 6 

Same  as  Fig.  5  for  standard  CaS  phosphor.  Compare  with  Figs.  2  and  4 
Fig.  7 

Spectral  distribution  of  IR  photons  capable  of  stimulating  the  standard  MgS 
phosphor  after  UV  excitation.  The  intensity  is  a  relative  measure  of  the 
effectiveness  of  a  given  wavelength  for  IR  stimulation. 
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studied  the  spectal  dependence  of  the  decay  after  UV  excitation  with  the  time 
resolution  of  approximately  300  ms  possible  with  the  OMA  III  system  in  its 
present  configuration.  Figure  9  indicates  that  the  Ce  luminescence  has 
completely  decayed  in  less  than  300  ms.  A  gated  OMA  detector  must  be  employed 
to  improve  the  time  resolution.  This  option  is  presently  not  available  in  our 
laboratory. 

For  applications  where  weaker  excitation  sources  must  be  employed  or  storage 
of  the  latent  image  is  desired,  the  storage  time  after  excitation  must  be 
known.  We  have  studied  the  response  to  a  given  IR  stimulation  radiation  after 
excitation  for  a  fixed  time  with  a  constant  UV  source.  The  results  shown  in 
Figs.  10  and  11  indicate  that  at  room  temperature  storage  times  up  to  one  hour 
(and  possibly  much  longer)  are  possible  without  significant  loss  of  the  stored 
image  information. 
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Fig.3 

Same  as  Fig. 7  but  for  standard  CaS  phosphor 

Fig.  9 

Decay  of  the  luminescence  emission  from  the  standard  MgS  phosphor  after  UV 
excitation.  The  Ce  emission  has  completely  decayed  after  less  than  300  ms. 

Fig.  10 

Test  of  storagebility  at  room  temperature  of  the  energy  deposited  during  UV 
excitation.  The  standard  MgS  phosphor  was  excited  for  a  fixed  time  with  253.7 
nm  UV  photons.  The  luminescence  emission  during  IR  stimulation  was  measured 
after  the  storage  time  indicated  on  the  curves,  which  are  offset  for  easier 
identification.  The  sharp  rise  of  these  response  curves  coincides  with  the 
onset  of  exposure  to  the  broadband  IR  stimulation  radiation.  The  increase  in 
signal  up  to  30  min  storage  time  is  real  and  reproducible,  kt  room  temperature 
the  stored  signal  decreases  only  after  storage  times  of  60  min.  or  more. 

Fig.  11 

Same  as  in  Fig.  10  but  for  the  standard  CaS  phosphor.  'lote  the  significant 
decrease  of  the  stored  signal  after  only  a  few  minutes  of  storage  at  room 
temperature. 
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